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5-H ■ ABSTRACT 

ryj I We study the 37 brightest radio sources in the Subaru/ XMM-Newton Deep 

. Field (SXDF). We have spectroscopic redshifts for 24 of 37 objects and 

photometric redshifts for the remainder, yielding a median redshift z me d 
for the whole sample of z me d — 1.1 and a median radio luminosity close 
to the 'FRI/FRII' luminosity divide. Using mid-IR (Spitzer MIPS 24 /mi) 
data we expect to trace nuclear accretion activity, even if it is obscured 
at optical wavelengths, unless the obscuring column is extreme. Our re- 
OO ■ suits suggest that above the FRI/FRII radio luminosity break most of the 

I radio sources are associated with objects that have excess mid-IR emis- 

, sion, only some of which are broad-line objects, although there is one clear 

C^) ' low-accretion-rate object with an FRI radio structure. For extended steep- 

OO ! spectrum radio sources, the fraction of objects with mid-IR excess drops 

■ dramatically below the FRI/FRII luminosity break, although there ex- 
ists at least one high-accretion-rate 'radio-quiet' QSO. We have therefore 
shown that the strong link between radio luminosity (or radio structure) 

■ and accretion properties, well known at z ~ 0.1, persists to z ~ 1. In- 
vestigation of mid-IR and blue excesses shows that they are correlated as 
predicted by a model in which, when significant accretion exists, a torus 
of dust absorbs ^30% of the light, and the dust above and below the torus 
scatters >1% of the light. 
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1 INTRODUCTION 

There is growing evidence that all massive galax- 
ies contain a supermassive black hole (SMBH) at 
their centres (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000; Magorrian 2006), which implies that these 
galaxies have all undergone at least one episode of 
AGN activity. It is believed that this AGN activity is 
the observational manifestation of feedback processes 
regulating galaxy and cluster evolution (Best et al. 
2005). At low redshifts radio jets emerge from the 



vicinities of these black holes, where there are indi- 
cations for a 'duty cycle' for active jet activity that 
varies from ~ 0.1 — > 0.3 over the 1.4-GHz luminos- 
ity range Li. 4 gh z ~ 8 x 10 23 ^8x 10 21 WHz-'sr" 1 
- a critical range as it provides a significant fraction 
(~25%) of the luminosity density of the population. 

Powerful (L1.4GHZ > 8 x 10 23 WHz _1 sr _1 ) radio 
sources are more common at high (z>l) redshifts. 
They are believed to have central SMBHs with uni- 
formly high accretion rates at the highest radio lu- 
minosities and relatively low duty cycles (e.g. Rawl- 
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ings & Saunders 1991). Low-luminosity radio jets can, 
however, be associated with high-accretion-rate sys- 
tems, and these so-called 'radio-quiet' quasars appear 
to have similar FRI-like radio structures to the low- 
accretion-rate counterparts of similar radio luminos- 
ity objects (e.g. Heywood et al. 2007). At low redshift, 
the most massive (> 1O S M ) SMBHs typically have 
very low accretion rates with systematically higher 
average values at z> 2, the so-called 'quasar epoch' 
(Yu & Tremaine 2002). These observational results 
fit in with theoretical ideas that a 'quasar mode' of 
feedback is prevalent in the distant universe, and that 
a 'radio mode' of feedback is dominant at low redshift 
(e.g. Croton et al. 2006). 

In the unified model for AGNs (e.g. Antonucci 
1993), the central region is surrounded by a dusty 
torus which absorbs light and re-emits it in the in- 
frared. Above and below the plane of the torus, 
dust scatters light yielding a blue excess (di Serego 
Alighieri et al. 1993; Tadhunter et al. 1992; Cimatti et 
al. 1993; Jannuzi et al. 1995). Such mechanisms make 
it difficult to observe objects viewed through the torus 
directly in the optical, UV and soft X-rays. The torus 
creates anisotropic obscuration of the central regions 
resulting in two different types of observed objects, 
type 1 that are viewed face-on and type 2 that are 
viewed edge-on. Here we use mid-IR observations to 
search for evidence of accretion in a manner which 
is far less dependent on orientation. Note, however, 
that there are claims that hot-dust emission is mildly 
anisotropic due to the effects of high optical depths 
and a toroidal geometry (Granato & Danese 1994). 

Through the years, there have been many at- 
tempts to classify radio sources according to their ob- 
served properties: e.g. their structure, energy output 
and environment. Fanaroff & Riley (1974) classified 
radio sources according to their structure by measur- 
ing the ratio of distance between the regions of high- 
est surface brightness on opposite sides of the central 
galaxy or quasar to the total extent of the source up 
to the lowest brightness contour Rfr- This way, for 
FRI radio sources we get -Rfr < 0.5 which makes 
them edge-darkened, while FRIIs have Rfr > 0.5 
meaning they are edge-brightened. Furthermore, they 
found a strong correlation of structure with radio 
luminosity: above log 10 (Li78MHz/WHz _1 sr _1 ) Ftf 25 
{H = 50 km s _1 Mpc _1 , fi M = 1 and Q A = 0)* lie 
the FRIIs, and below that the FRIs. 

Ledlow & Owen (1996) found that both FRIs and 
FRIIs live in similar environments, but the properties 
of their host galaxies may influence the structural ap- 
pearance of the radio sources, at least for those near 
the division of the two FR classes. A more refined 



* This corresponds to log 10 (Li7gMHz/WHz _1 sr~ 1 ) fid 25 
in our adopted cosmology (Hq = 70 km s _1 Mpc _1 , CIm = 
0.3 and Q A = 0.3) at the typical redshift (z ~ 0.1) of the 
FRI/FRII division from the Fanaroff & Riley (1974) study. 



division of radio sources by structure is the one ap- 
pearing in Owen & Laing (1989). In this case we have 
Classical Double (CD), Twin Jet (TJ) and Fat Dou- 
ble (FD) radio sources: CDs correspond to FRIIs, TJs 
to FRIs and FDs to FRI/FRII division objects. Study 
of the 3CRR, 6CE, 7CRS and TOOT radio samples 
at 2 ~ 0.5 (McLure et al. 2004) is consistent with 
the low-redshift result of Owen & Laing (1989) that 
FRI (FD & TJ) sources reside in hosts which are on 
average ~ 0.5 mag brighter than those of classical 
double/FRII sources of comparable radio luminosity. 

Other classification schemes, such as the 'dual- 
population' model (e.g. Jackson & Wall 2001), di- 
vide objects according to nuclear accretion rate fol- 
lowing the pioneering study of Hine & Longair 
(1979). The less radio-luminous population is com- 
posed of FRIs and FRIIs with weak or absent nar- 
row emission lines, and the more radio-luminous 
population of strong narrow-emission-line FRII ra- 
dio galaxies and broad-line quasars, where the divi- 
sion is at log 10 (Li7gMHz/WHz _1 sr _1 ) » 26 {Ho = 
70 km s^Mpc -1 , fi M = 0.3 and fi A = 0.7) and cor- 
responds to the Radio Luminosity Function (RLF) 
break (e.g. Dunlop & Peacock 1990). Best et al. 
(2003) investigated the dual-population scheme for 
radio sources using the CENSORS sample, selected 
from the 1.4 GHz NVSS survey (Condon et al. 1998). 

Ogle et al. (2006) studied a sample of FRII 
narrow-line radio galaxies from the 3CRR survey at 
redshifts z < 1, in order to investigate whether they 
host hidden quasar nuclei. They found that in most 
of these radio galaxies, absorption in the mid-IR is 
present due to dust from a molecular torus. Almost 
half of their sample of narrow-line FR II radio galaxies 
is weak at 15 or 30 /im, contrary to single-population 
unification schemes. They conclude that mid-IR weak 
radio galaxies may constitute a separate population 
of jet-dominated sources with low accretion power. 
In Vardoulaki et al. (2006) we expressed the need to 
study radio sources in the mid-IR in order to investi- 
gate the existence of any hidden accretion activity in 
low-luminosity radio sources. Spitzer observations at 
24 fim are ideal for this task, since they trace warm 
dust emission which can be obscured only by extreme 
columns (Ogle et al. 2006). 

To address this important issue, we have looked 
at a complete sample of radio sources in the 
Subaru/ XMM-Newton Deep Field (SXDF). The sam- 
ple studied here is the 37 brightest radio sources from 
the VLA survey (SXDS) of the SXDF (Simpson et 
al. 2006) with flux densities greater than 2 mjy at 
1.4 GHz. In Section 2 we present new data taken on 
the SXDS radio sources, namely optical and near- 
infrared imaging and spectroscopy. We also present 
a cross-correlation with the radio sample of Tasse et 
al. (2006), who observed this sky region at 325 MHz, 
as well as GMRT observations at 610 MHz (Sec. 2.1). 
Section 3 presents the analysis of our sample including 
if-band/radio overlays (Sec. 3.1), photometric red- 
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shift estimation (Sec. 3.2), and the classification of the 
SXDS-based sample according to optical and infrared 
data (Sec. 3.3). In Section 4 we discuss accretion in- 
dicators, based on our optical/infrared radio-source 
classification. In Section 5 we present the conclusions 
of our study. Finally, in the Appendix, we describe 
properties of each radio source individually. 

We adopt a radio spectral index a — 0.8 (S„ oc 
v~ a ), unless a spectral index could be calculated us- 
ing the 1.4 GHz, 610 MHz and 325 MHz data. We 
assume throughout a low-density, A-dominated Uni- 
verse in which Ho = 70 km s _1 Mpc _1 , Qm = 0.3 and 
Ov = 0.7. 



2 DATA AND OBSERVATIONS 

Optical and X-ray observations of the SXDF 
(Fig. 1) were made within the 1.3 square degree 
Subaru/ XMM-Newton Deep Field with Subaru and 
XMM-Newton respectively. A follow-up deep radio 
imaging survey at 1.4 GHz was conducted with the 
Very Large Array (VLA) in B and C configurations 
(Simpson et al. 2006). A-array observations are also 
available (Ivison et al. 2007). The total radio cata- 
logue, hereafter the SXDS, covers 0.8 deg 2 to a peak 
flux density limit of 100 /iJy at 1.4 GHz. 

The Subaru-XMM-Newton Deep Field (SXDF) 
has excellent multi-wavelength coverage, having been 
observed in X-rays with the XMM-Newton satel- 
lite (Ueda et al., in prep.), in the optical BVRi'z' 
bands with Suprime-Cam on the Subaru telescope 
(Sekiguchi 2004), in the near infrared J and K bands 
as part of the UKIRT Infrared Deep Sky Survey 
Ultra-Deep Survey (UKIDSS - UDS; Lawrence et al. 
2007), in all of the Spitzer IRAC (3.6, 4.5, 5.8 and 
8.0 fj,m) and MIPS (24, 70 and 160 fj,m) bands as 
part of the Spitzer Wide-area Infra-Red Extragalactic 
(SWIRE) survey Lonsdale et al. (2003) and at sub- 
mm wavelengths (450 /im and 850 /im; Mortier et al. 
2006 & Aretxaga et al. 2007). The area is large (for 
such a deep survey) , at approximately 1 square degree 
in the optical and near-infrared, although the cover- 
age varies between the different wavebands. The UDS 
covers approximately 0.8 deg 2 centred on the SXDF. 

We study here the brightest 37 radio sources from 
the 1.4 GHz SXDS (Simpson et al. 2006); these are 
named, hereafter, sxds_0001 to sxds_0037 in decreas- 
ing order of 1.4-GHz flux density. Table 1 gives basic 
information about our sample, like their radio posi- 
tion, flux density, redshift, optical and radio classifi- 
cation, spectral index, angular size and K (2.2 /im) 
magnitude (in the AB magnitude system). We adopt 
the radio structure classification from Owen & Laing 
(1989), using the VLA data to classify objects as 
CD, TJ, FD or COM (compact) radio sources. Wher- 
ever available, we used the A-array data to determine 
whether an object has compact hotspots and is there- 
fore a 'CD': in the few cases where this classification 



is uncertain, because of no available A-array data, we 
took our best guess based on the appearance of the 
B-array map and marked the classification with a '?'. 
Details of the optical classification are presented in 
Section 3.3. 

In Table 2 we give the result of optical spec- 
troscopy of the SXDS radio sources. We have w 65% 
(24 out of 37) completeness in spectroscopic red- 
shifts, with the spectra shown in Fig. Al. Table 2 
summarises spectroscopic observations carried out on 
the SXDS objects using either the ISIS spectrograph 
on the WHT telescope, FOCAS on the SUBARU, 
2dF or AAOmega on the AAT, VIMOS on the ESO- 
VLT, SITel on the Magellan 2 or DEIMOS on the 
Keck 2. Calculation of redshifts from VIMOS spec- 
tra were done by using the fxcor cross-correlation 
task in IRAF with SDSS galaxy templates (Simp- 
son et al. in prep), apart from sxds_0025 where the 
cross-correlation failed (see Appendix A). FOCAS 
spectra have been corrected for atmospheric absorp- 
tion. The object sxds_0016 has spectroscopy obtained 
with the DEep Imaging Multi-Object Spectrograph 
(DEIMOS) on the Keck 2 telescope (van Breukelen 
et al. 2007). 

Table 3 presents the blueness and mid-IR excess, 
the radio luminosity at 1.4 GHz Li.4ghz, the lumi- 
nosity at 24 /im [AL]24 M m and basic quantities for the 
SXDS radio sources used in the analysis (see Sec. 4.1). 

Optical and infrared photometric data were used 
to create SEDs for our objects as shown in Fig. Al. 
Thirteen of our objects are not as yet spectroscopi- 
cally confirmed, so we use photometric redshifts for 
our analysis (Sec. 3.2). In Figure 1 we have shown 
the sky area occupied by the 37 radio sources in our 
sample. We are confident that our sample is fairly 
complete. A NVSS-selected 2-mJy sample in the same 
area gets 34 of the 37 SXDS sources, and the other 
four lie just below the formal NVSS 2.3 mjy flux den- 
sity limit. 

In order to calculate spectral indices for our sam- 
ple we use either the 325 MHz catalogue of Tasse et 
al. (2006), or the 610 MHz flux densities from Ibar 
et al. (in prep), taking into account the radio-source 
structure from the A/B-array maps. 

Tassc et al. (2006) surveyed the XMM-LSS field 
(Pierre et al. 2004) in A- and B-array VLA configura- 
tions at 325 MHz (see Sec. 6). To minimise the effect 
of Doppler-boosted emission, we would ideally have 
selected our objects at 325 MHz, but this survey was 
unavailable at the time the selection was made. We 
have, however, subsequently cross-matched the 325 
MHz sample with the whole 1.4 GHz radio catalogue. 
The resolution of the 325 MHz data is 6.7 arcsec and 
it has a median 5a sensitivity of 4 mjy beam -1 . Fig- 
ure 1 shows the 325-MHz sources that are in a sky 
area somewhat larger than our SXDS sample. From 
our 37-source sample only 22 match with the 325 MHz 
catalogue. For all of the matching objects we calculate 
radio spectral indices using the 1.4 GHz and 325 MHz 
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Figure 1. Sky positions of the 37 brightest radio sources in the SXDS (blue circles); these are scaled by the logarithm of 
the integrated flux density at 1.4 GHz from the VLA observations of Simpson et al. (2006) starting at Si.4GHz = 2 mjy. 
The shaded area shows the UKIDSS-UDS (Lawrence et al. 2007) coverage. Note that sxds_0009, sxds_0017, sxds_0025, 
sxds_0027 and sxds_0037 lie outside the UDS coverage. Light-blue stars are the VLA A-array pointing positions and big 
light-blue circles the FWHM points of the primary beam of the VLA primary beam; objects inside these regions have 
reasonable quality A-array maps and are highlighted with filled blue circles. Green triangles are the 325 MHz sources 
(S l 325MHz~ 4 mjy) in this area (Tasse et al. 2006). Red crosses are NVSS (Condon et al. 1998) sources in the same sky 
area, again with the symbol size scaled by the logarithm of the NVSS flux density starting at 2 mjy. The 1 Mpc bar shows 
the proper physical size of a 1 Mpc-long ruler at z = 1. 



flux densities. Limits on the spectral indices are cal- 
culated assuming a S325MHz~ 4 mjy (5a) limit for the 
325 MHz catalogue, although (see Appendix A) there 
are a few cases where this process is unreliable (e.g. 
sxds_0011 where a steep spectral index is expected 
from the TJ radio structure of this radio source). 

The fact that there are 13 325-MHz objects with 
reasonably high flux density that match with lower 
1.4-GHz flux density objects than our 2 mjy limit 
is unsurprising, given that S325MHz~ 4 mjy corre- 
sponds to Si.4GHzS; 2 mjy for a ~ 0.9 and such steep- 
spectrum objects exist; their number will be exagger- 
ated by random biases on the 325-MHz flux densities 
at the limit of the Tasse et al. (2006) survey. 

The 610-MHz observations of the SXDF were 
made with the Giant Metre-wave Radio Telescope 



(GMRT) in two observing runs: on 2006 February 2-5 
and 2006 December 4-9 (details can be found in Ibar 
et al. in prep). The region was covered with three 
pointings forming an equilateral triangle with 22- 
arcmin sides, covering an area of ~0.5 deg 2 . Adding 
both sets of observations, the total integration time 
per pointing was 11.4 hr. 

To calibrate the data, we used 3C 147 or 3C 48 
as flux and bandpass calibrators, observed for 20 min 
at the beginning and end of each track. We also ob- 
served 0116-208 or 0240-231 for 7 min every 45 min 
to be used as phase and amplitude calibrators. The 
full available bandwidth of 32 MHz was split into two 
intermediate frequencies (IFs) of 128 channels each, 
and reduced independently. 

For the data reduction we used the Astronomi- 
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cal Image Processing System (AIPS). Standard tasks 
(e.g. QUACK, tvflg, UVFLG, SPFLG) were used to 
flag bad baselines, antennas, initial integrations, and 
channels that suffered from narrow-band interference. 
A first rough calibration using CALIB and CLCAL was 
performed before using the flgit task to flag more 
bad data using a noise-based criterion. The clean data 
were calibrated again and used to normalise the tar- 
get using GETJY. Each IF was compressed in chunks 
of 7 channels (to ensure bandwidth smearing was not 
a problem) leaving 15 averaged channels which were 
used for imaging. Each pointing was broken into 37 
facets, in order to minimise 3D non-coplanar issues, 
and imaged using IMAGR. 

The GMRT has a large number of antennas in 
the central region. These dominate the uv coverage 
and result in a non-Gaussian beam shape. Baselines 
shorter than 1.5 kA were omitted from the imaging. 
We utilised one self-calibration in phase and two in 
phase and amplitude, averaging both polarisations 
in these two last passes, in order to obtain the fi- 
nal image facets. These were combined using flatn, 
after convolving them to a single, common beam size. 
The final radio image reaches an r.m.s. of ~ 60 /iJy 
beam" 1 near the centre of the field, with a synthesised 
beam of 6.8 x 5.4 arcsec 2 at position angle of 30 deg. 
The south-west field, centred at R.A. 2h 18m 44.0s, 
Dec. —5° 7' 20" has a considerably larger mean r.m.s. 
(~ lOO^iJy) than the other two fields. Also, the noise 
in the image increases near bright sources, where it 
can be as high as ~ 300 /iJy r.m.s., limiting the de- 
tection of nearby sources. 



3 ANALYSIS 

3.1 Basic information 

In our analysis we use the radio flux densities from Ta- 
ble 1 according to the following criteria. In the cases 
where the NVSS flux density is higher than the SXDS 
one, and the source is not compact, we use the NVSS 
flux density for our analysis, since the source may 
have more structure than the B-configuration VLA 
observations could reveal. On the other hand, when 
the NVSS flux density is higher than the SXDS flux 
density and the source is compact we assume that the 
NVSS measurement is inaccurate due to confusion by 
a nearby source (which we confirmed by inspection of 
the VLA image); in this case we use the SXDS flux 
density. Finally, when the SXDS flux density is higher 
than the NVSS value, and the source has compact ra- 
dio structure, we assume that the source is time vari- 
able; in that case we also use the SXDS flux density 
for our analysis. 

For every SXDS radio source if— band UDS im- 
ages were downloaded from the WFCAM archive 
(Hambly et al. 2007), with the exception of sxds_0009, 
sxds_0017, sxds_ 0025, sxds_0027 and sxds_0037 that 



Table 1. Basic properties of the 37 brightest SXDS ra- 
dio sources; values in bold are those used in the analy- 
sis. Columns 1, 2 & 3 give the name of the object and 
its VLA (1.4 GHz) radio position (J2000.0) from Simp- 
son et al. (2006). Columns 4 & 5 give the integrated 
flux densities at 1.4 GHz from Simpson et al. (2006) and 
NVSS respectively; radio sources labelled with 'x' suffer 
from confusion by a nearby source, which results in a 
higher NVSS flux density than the SXDS value as seen 
in objects sxds_0025 and sxds_0027. Objects sxds_0034, 
sxds_0036 and sxds_0037 are not in the NVSS catalogue 
since they are fainter than the 2.3 mjy limit of the survey 
(Condon et al. 1998); their flux densities were measured 
from the NVSS maps and are labelled with 'n'. Column 6 
gives the 610 MHz flux density in mjy (Ibar et al. in prep); 
the character 's' shows that the flux densities of each com- 
ponent of the radio source have been added to calculate 
the whole flux density of the radio sources; '(nif)' means 
that the object is not inside the field that was observed; 
values in brackets are inaccurate since they miss part of 
flux density coming from the source. In Column 7 we 
give the 325-MHz flux density from Tasse et al. (2006); 
values in brackets are not accurate since they miss part 
of the flux density coming from the source. Column 8 
shows our radio source counterpart optical/IR classifica- 
tion (see Sec. 3.3): 'Q' stands for a quasar, 'OQ' for an 
obscured quasar, 'G?' for a possible galaxy, 'G' for a se- 
cure galaxy, 'SB' for a starburst, 'WQ' for a weak quasar 
and 'BL' for a BL Lac object. Column 9 gives the spec- 
troscopic redshift, when available: '?' denotes objects with 
spectroscopy yielding plausible but insecure rcdshifts; the 
symbol '+' indicates clear broad emission lines and 'a' de- 
notes a broad-absorption-linc QSO (BALQSO). Columns 
10 & 11 give the photometric redshift calculated with HY- 
PERz using Coleman (CWW) and Bruzual-Charlot (BC) 
templates respectively; brackets are used to denote unreli- 
able fits. Column 12 gives the radio spectral index from 
1.4 GHz to 610 MHz or from 1.4 GHz to 325 MHz, which 
are marked with 'm'; we use the 1.4 GHz flux densities in 
bold for the calculations. Estimated values are in brack- 
ets. Column 13 gives the largest angular size in units of 
arcsec and Column 14 gives the radio classification (R 
CI) of each object (Owen & Laing 1989): CD stands for 
classical double; TJ for twin-jet; FD for fat double; and 
COM for a compact radio source. 



were not observed by the UDS (Fig. 1). For these ob- 
jects we use z— band images from Subaru in order to 
produce overlays of the radio B-array VLA data and 
imaging, as shown in Fig. Al. Also in Fig. Al we 
present spectral energy distributions (SEDs) of the 
objects created as explained in Section 3.2. 
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02 


16 


16 


.82 


-05 


12 


53 


.47 


6.25 


6.9 


6.7 


5.4 


sxds_0014 


02 


17 


52 


.14 


-05 


05 


21 


.25 


6.19 


5.7 


4.2 




sxds_0015 


02 


19 


32 


.20 


-05 


07 


32 


.66 


6.10 


6.0 


17.0 s 




sxds_0016 


02 


18 


26 


.12 


-04 


59 


46 


.57 


5.95 


4.4 


14.1 s 


16.9 


sxds_0017 


02 


18 


27 


.57 


-05 


34 


53 


.77 


4.97 


4.6 


(nif) 


13.5 


sxds_0018 


02 


17 


24 


.38 


-05 


12 


51 


.68 


4.84 


5.7 


13.3 




sxds_0019 


02 


17 


57 


.26 


-05 


27 


55 


.82 


4.83 


5.0 


(nif) 




sxds_0020 


02 


18 


00 


.68 


-04 


49 


54 


.78 


4.60 


5.1 


8.5 s 


3.9 


sxds_0021 


02 


17 


52 


.53 


-04 


48 


23 


.90 


4.42 


5.0 X 


9.1 


10.2 


sxds_0022 


02 


17 


55 


.15 


-05 


26 


52 


.99 


4.36 


4.7 X 


6.1 


9.5 


sxds_0023 


02 


17 


54 


.10 


-05 


12 


49 


.94 


4.20 


4.1 


11.3 


15.2 


sxds_0024 


02 


19 


06 


.64 


-04 


59 


00 


.56 


3.61 


3.8 


11.2 


4.02 


sxds_0025 


02 


17 


55 


.42 


-05 


37 


04 


.97 


3.23 


5.9 X 


(nif) 




sxds_0026 


02 


18 


57 


.04 


-05 


28 


19 


.98 


3.20 


3.4 


6.4 




sxds_0027 


02 


17 


54 


.93 


-05 


36 


29 


.10 


3.13 


5.9 X 


(nif) 




sxds_0028 


02 


17 


18 


.59 


-05 


29 


21 


.02 


2.80 


5.2 X 


4.1 


_ 


sxds_0029 


02 


18 


39 


.58 


-05 


01 


34 


.68 


2.72 


4.4 X 


3.4 




sxds_0030 


02 


18 


22 


.44 


-05 


16 


48 


.77 


2.56 


2.4 


4.6 


5.4 


sxds_0031 


02 


18 


22 


.65 


-05 


02 


53 


.26 


2.44 


3.2 


3.9 




sxds_0032 


02 


19 


26 


.48 


-05 


15 


35 


.00 


2.39 


2.3 


7.7 s 




sxds_0033 


02 


17 


37 


.21 


-05 


13 


27 


.96 


2.37 


2.4 


6.3 


3.9 


sxds_0034 


02 


18 


09 


.48 


-04 


59 


45 


.54 


2.31 


1.21™ 


4.2 


6.2 


sxds_0035 


02 


16 


59 


.52 


-05 


16 


55 


.73 


2.24 


2.9 X 


1.5 




sxds_0036 


02 


19 


32 


.15 


-05 


12 


57 


.26 


2.09 


2.72™ 


<1.1 




sxds_0037 


02 


18 


06 


.88 


-05 


37 


17 


.67 


2.07 


1.68™ 


(nif) 





(8) (9) (10) (11) (12) (13) (14) 



Dptical/IR 


Zspec 


z phot 


z phot 


a 


9 


R CI 


Type 




CWW 


BC 




(arcsec) 




G 


0.626 


0.618 


0.710 


1.03 


30 


TJ 


OQ 


1.183 


1.080 


1.020 


0.93 


4.4 


CD 


OQ 


2.433 


[2.0441 

L J 


[1.4451 

L J 


1.17 


9 


CD 


BL 


? 


[2.126] 


2.270 


0.34 


<2 


COM 


G? 


? 


1.633 


1.610 


1.07 


<2 


COM 


G? 




1.650 


1.882 


-0.65 


<2 


COM 


G? 


1.325? 


[1.174] 


0.980 


1.81 


123 


CD 


SB 


0.458 


[0.218] 


0.500 


0.72 


<2 


COM 


Q? 


3.570? 


[3.565] 


3.585 


0.53 


2.5 


COM 


G? 




[1.231] 


1.400 


0.90 


3.3 


COM? 


G 


0.647 


0.618 


0.695 


(0.8) 


120 


TJ 


G 


0.865 


0.868 


0.770 


0.97 


120 


CD 


Q 


2.713+ a 


[2.922] 


[2.475] 


0.08 


<2 


COM 


G? 


? 


[0.240] 


2.205 


-0.47 


<2 


COM 


G 


0.344 


0.386 


0.385 


1.23 


47 


FD 


OQ 


1.345? 


1.008 


1.082 


1.04 


45 


CD 


Q 


2.583+ a 


[2.771] 


[2.120] 


0.68 m 


22 


CD? 


G 


0.919 


0.858 


0.820 


1.02 


100 


FD 


G 


0.694 


[2.071] 


0.910 


(0.8) 


52 


FD? 


G? 




[2.553] 


1.135 


0.62 


35 


FD? 


G? 




1.728 


1.645 


0.87 


<2 


COM 


G 


0.248 


0.280 


0.290 


0.40 


<2 


COM 


G 


0.586 


0.561 


0.610 


1.19 


3 


COM 


G 


0.516 


0.459 


0.500 


1.30 


45 


TJ 


WQ? 


0.964? 


[0.019] 


[2.465] 


(0.4) 


<2 


COM 


G 


0.450 


0.430 


0.445 


0.76 


10 


TJ? 


OQ? 


? 


[0.011] 


1.595 


(0.4) 


<2 


COM 


G 


0.627? 


[1.092] 


0.895 


0.46 


<2 


COM 


G? 




1.272 


1.380 


0.27 


<2 


COM 


G 


0.535 


0.511 


0.650 


0.71 


4 


TJ? 


G? 




1.335 


1.325 


0.24 


<2 


COM 


G 




1.196 


1.010 


1.41 


15 


FD? 


G 


0.647 


0.774 


0.790 


1.16 


18 


TJ 


OQ 


1.095 


[1.309] 


[1.345] 


0.72 


<2 


COM? 


G? 




3.144 


2.715 


-0.48 


<2 


COM 


G 


0.871 


0.857 


0.840 


<-0.77 


<2 


COM 


G? 




[0.709] 


1.085 


(0.4) 


<2 


COM 
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Table 2. Observing log for optical spectroscopy of the SXDS radio sources: Column 1 gives the name of the object. 
Column 2 gives the optical position of the radio sources in RA (h m s) and DEC (° m s). Column 3 shows the instrument 
and telescope where spectroscopic observations were conducted; cither with ISIS at the WHT (www.ing.ias.es), with 
FOCAS at the SUBARU (www .subarutelescope.org) , with VIMOS at the ESO-VLT (www .eso.org / instruments /vimos /) , 
with 2dF/AAOmega at the AAT (www.aao.gov.au/AAO/), with SITel at the Magellan 2 (www.ociw.edu/Magellan), or 
with DEIMOS at Keck 2 (www.keckobservatory.org). Column 4 gives the observing date and Column 5 the spectroscopic 
rcdshift (a '?' is added when the redshift not reliable). Column 6 gives the exposure time; in the case of the WHT we 
give exposure time for both the ('B') and ('R') parts of the spectrum. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Object 


Optical position 


Spectrometer/Telescope 


Date 


^spec 


Exposure 












Time (s) 


sxds_UUUl 


no 1 Q 07 1 ftl 

yjz lo z i . ibi 


V IIVIUo/HjoU- V Lil 


ocnp+ onnfc 
zbuctzuuo 


O.DZD 


o4UU 




-04 54 41.60 


TflTfl /T T TT TJT1 

ISIS/ WHT 


26Jul2003 


0.625 


lx900B, 2x900R 


sxdsJJUUz 


no 1 Q IS 1 'JO 

VZ lo lo.loy 


lolo- Wli 1 


of; TnionnQ 


1 1 QQ 


IxloUUrj, zxyUUrt 




f\ A AC r\*7 EC EC 

-U4 4b U7.55 










sxds_0003 


02 18 39.559 


VIMOb/bbO- VL1 


OOZ~\ j-OOOE" 

28Oct2005 


2.432 


5400 




C\A /II /in X.7 
-U4 41 4y.O ( 


lolo/ V V XT 1 


os TnionnQ 




IxloUUrj, zxyUUrt 


sxas_uUU4 


no 1 Q EC Q EC Q "7 

U2 lo 53. oo f 


A A Pil T T? r" 1 A / A AT 1 

A AU M CjKj A/ A A 1 






1800 




-04 47 36.29 










sxas_UUUo 


no 1 q cc i qoc 
U2 lo 51.32b 


TQTQ /WJUT 
lolo/ V V XT 1 


07 t,, i qooq 
I ( JulzUU3 


9 


IxloOUrj, zxyOUrt 


sxds_0006 


02 lb 37.825 












ncc. 1 EC oa os 

-UO lO zo.ZO 










sxds_0007 


02 lb 59.064 


1 OOAb/bUBARU 


22Dec2003 


1.325 


1200 




H/1 /in on s ec 
-U4 49 zU.oO 










sxds_0008 


oo 1 o o .i nnr 

02 18 24.005 


O J TT* / A A rn 

2dl /AA1 


7Dec2002 


0.457 


1800 




-04 53 05.23 


VIMOb/EbO-vTT 


26Oct2005 


0.458 


5400 


sxds_0009 


02 18 03.420 


A A /~~\ A /rT7l/"1 A / A A rri 

AAOMbGA/AAl 




3.570? 


1800 




of O O Of JO 

-05 38 25.43 










sxds_0010 


02 18 50.491 












n/1 ECO QO Q1 

-U4 Do oz.ol 










sxds.OOll 


o o 10 oo con 

02 18 23.532 


O J T7I / A A Ti 

2dl /AA1 


7Dec2002 


0.647 


1800 




aec o ec oo 

-05 25 00.69 


bl lei /Magellan 2 


18Uct2003 


0.647 


2280 


sxds_0012 


02 16 34.968 


VIMOb/EbO-VLT 


8Aug2005 


0.865 


5400 




OA flfi 47 
-U4 OO UD.4f 










sxds_0013 


02 16 16.822 


FOCAb/bUBARU 


6Oct2003 


2.713 


1800 




-05 12 53.71 










sxds.0014 


02 17 52.118 


AAOMEGA/AAT 




? 


1800 




-05 05 22.23 










sxds.0015 


02 19 32.201 


2dF / AAT 


12Dec2002 


0.344 


1800 




-05 07 32.67 


AAOMEGA/AAT 




0.344 


1800 


sxds.0016 


02 18 26.116 


DEIMOS/KECK 2 


21Dec2006 


1.345? 






-04 59 46.57 










sxds.0017 


02 18 27.312 


FOCAb/bUBARU 


23Dec2003 


2.583 


1800 




-05 34 57.40 










sxds.0018 


02 17 24.415 


FOCAb/bUBARU 


28Oct2003 


0.919 


1800 




-05 12 52.63 










sxds.0019 


02 17 57.288 


FOCAb/bUBARU 


20Dec2003 


0.694 


1800 




-05 27 55.88 










sxds_0020 


02 18 00.716 












-04 49 56.41 










sxds.0021 


02 17 52.531 












-04 48 24.68 










sxds.0022 


02 18 18.139 


VIMOb/EbO-VLT 


2bep2005 


0.248 


5400 




-04 46 07.55 










sxds.0023 


02 17 54.125 


VIMOS/EbO-VLT 


8Oct2005 


0.586 


5400 




-05 12 50.14 










sxds.0024 


02 19 06.665 


VIMOb/EbO-VLT 


10Oct2005 


0.516 


5400 



-04 59 01.78 
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Table 2. (continued) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Object 


Optical position 


Spectrometer/Telescope 


Date 


■^spec 


Exposure 
Time (s) 


sxds_0025 


02 17 55.392 
-05 37 05.00 


VIMOS/ESO-VLT 


7Oct2005 


0.936? 


5400 


sxds_0026 


02 18 57.046 


2dF / AAT 


8Dec2002 


0.450 


1200 




-05 28 19.38 


2dF / AAT 


10Dec2002 


0.450 


1200 


sxds_0027 


02 17 54.943 
-05 36 29.80 


VIMOS/ESO-VLT 


7Oct2005 


? 


5400 


sxds_0028 


02 17 18.564 
-05 29 20.71 


FOCAS/Subaru 


25Dec2003 


0.627? 


1200 


sxds_0029 


02 18 39.586 
-05 01 35.63 










sxris 0030 


02 18 22.493 
-05 16 48.83 


AAOMEGA/AAT 


31 Aue-2006 


0.535 


1800 


sxds_0031 


02 18 22.620 
-05 02 53.77 










sxds_0032 


02 19 26.484 
-05 15 34.62 










sxds_0033 


02 17 37.198 


2dF/AAT 


10Dec2002 


0.647 


1800 




-05 13 29.58 


SITel/Magellan 2 


HMar2005 


0.647 


2400 


sxds_0034 


02 18 09.456 
-04 59 45.89 


FOCAS/SUBARU 


18Dec2003 


1.095 


1200 


sxds_0035 


02 16 59.546 
-05 16 56.32 










sxds_0036 


02 19 32.143 
-05 12 57.64 


FOCAS/Subaru 


20Dec2003 


0.871 


1200 


sxds_0037 


02 18 06.893 











-05 37 17.88 
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Table 3. Column 1 gives the name of the object. Columns 2 & 3 arc the blue excess (log in ([i^Ll ° /fi^liumrest)) 

u J 4U00Arest 

and the mid-IR excess (log w ([^L]io M mrest/[i / -£<]i M mrcst)) respectively. Columns 4 & 5 give the logarithms of the rest 
frame luminosity at 1.4 GHz and at 24 (im respectively. Column 6 is the Kab magnitude from the UDS; the character 
l K' denotes point sources at the if— band, whereas 'o' denotes point sources in the optical. For sxds_0009, sxds_0017, 
sxds_0025, sxds_0027 and sxds_0037 we give an estimate of the K magnitude from the SED of the object (see Fig.l); these 
values arc in brackets. 



v-v 


(2) 


(3) 


(4) 


(5) 


(6) 


Object 


Blueness 


mid-IR 


!ogio( 


l°gio( 


Kab 






excess 


il.4GHz/ 


[AL] 24/jm / 










[WHz-isr" 1 ]) 


[W]) 


(mag) 


sxds_0001 


-1.17 


< -1.40 


25.04 


< 36.88 


18.49 


sxds.0002 


-0.60 


-0.71 


25.57 


37.30 


19.93 


sxds_0003 


0.33 


< 0.58 


26.39 


37.98 


20.89 


sxds_0004 


-0.17 


0.88 


25.38 


38.04 


21.40^ 


sxds.0005 


-0.32 


< 0.04 


25.36 


< 37.72 


20.57 


sxds_0006 


-1.23 


< 0.17 


24.62 


< 38.16 


20.21 


sxds_0007 


-0.60 


< -0.58 


25.20 


< 37.64 


19.68 


sxds_0008 


-0.21 


-0.13 


23.71 


37.34 


18.24 


sxds_0009 


-0.96 


< 0.19 


25.61 


< 38.51 


(22.52)° 


sxds_0010 


-0.76 


< -0.39 


24.85 


< 37.60 


21.11 


sxds_0011 


-0.89 


< -1.51 


< 24.15 


< 36.89 


18.07 


sxds.0012 


-1.01 


< -1.31 


24.29 


< 37.08 


19.11 


sxds_0013 


-0.36 


< 0.73 


24.96 


< 38.25 


21.79° 


sxds_0014 


-3.19 


< -2.25 


24.51 


< 38.33 


21.23^ 


sxds.0015 


-0.33 


-0.85 


23.31 


36.66 


17.69 


sxds_0016 


-0.49 


< -0.19 


24.73 


37.86 


22.06 


sxds_0017 


0.30 


< 1.07 


25.15 


< 38.15 


(21.64)° 


sxds.0018 


-0.78 


< -1.22 


24.29 


< 37.18 


18.88 


sxds_0019 


-0.54 


< -0.57 


< 23.88 


< 36.97 


20.26 


sxds_0020 


0.06 


0.60 


24.33 


< 37.43 


21.22 


sxds.0021 


-0.09 


< 0.54 


24.80 


< 37.72 


21.04 


sxds_0022 


-0.54 


-1.54 


22.76 


35.87 


16.48 


sxds_0023 


-0.42 


< -1.02 


23.71 


< 36.71 


18.97 


sxds_0024 


-0.82 


< -1.63 


23.55 


< 36.66 


17.80 


sxds_0025 


0.09 


< 0.57 


< 23.91 


< 37.29 


(21.60) 


sxds_0026 


-1.00 


< -1.74 


23.27 


< 36.54 


17.65 


sxds_0027 


-0.39 


0.81 


< 24.37 


38.38 


(20.72) 


sxds_0028 


-0.27 


< -0.69 


23.45 


< 36.84 


18.71 


sxds_0029 


-0.91 


< -0.65 


24.12 


< 37.63 


20.62 


sxds_0030 


-0.85 


< -1.23 


23.31 


< 36.65 


18.44 


sxds_0031 


-0.78 


< -0.21 


24.15 


< 37.58 


20.32 


sxds.0032 


-0.47 


< -0.22 


24.13 


< 37.27 


20.51 


sxds_0033 


-0.23 


-1.02 


23.57 


36.83 


18.40 


sxds_0034 


-0.08 


1.00 


24.00 


38.60 


19.75^° 


sxds.0035 


-0.80 


< 1.50 


24.19 


< 38.28 


21.86 


sxds_0036 


-1.36 


< -1.32 


< 23.31 


< 37.08 


18.71 


sxds_0037 


-0.92 


< -0.52 


< 23.83 


< 37.40 


(20.60) 
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Table 4. Approximate 5-ct point-source depths for the 
bands used in the SXDF, and, for the Spitzer bands, the 
flux density limits of the SWIRE catalogues (used as upper 
limits in the SEDs for undetected objects). 



Band Depth Catalogue flux limit 





AB mags 








B 


27.5 


3.6 x 10- 


-2 




V 


26.7 


7.6 x 10" 


-2 




R 


26.8 


6.9 x 10" 


-2 




i' 


26.5 


9.1 x 10" 


-2 




z' 


25.5 


0.23 






J 


23.5 


1.4 






K 


23.5 


1.4 






3.6 fim 


22.5 


3.7 




10 


4.5 (an 


22.1 


5.4 




10 


5.8 fim 


19.7 


48 




43 


8.0 jitm 


20.0 


37.8 




40 


24 /im 


18.0 


230 




450 



3.2 Photometric redshifts 

We have estimated redshifts for the objects by fit- 
ting template spectral energy distributions - SEDs to 
their observed photometry, using a modified version of 
the publicly-available^ HYPERz code, by Bolzonella 
et al. (2000). HYPERz offers a convenient interface 
for the provision of photometric bandpass functions, 
galaxy templates and data, and is well documented. 
It also has built-in prescriptions for reddening of tem- 
plates by dust, and for Lyman forest absorption below 
912 A. 

For the purposes of galaxy template fitting, 
we consider only the optical (Subaru), near-infrared 
(UKIDSS-UDS) and mid-infrared (SWIRE) data. 
The approximate 5-cr point source depths achieved 
in each waveband are shown in Table 4. 

The catalogues contain a variety of adaptive and 
fixed aperture magnitudes, which have been mea- 
sured from the image data using SExtractor (Bertin 
& Arnouts 1996). For the purpose of SED fitting, it 
was important to select those apertures which pro- 
vided the most accurate colour information for each 
object. Whilst the adaptive (such as Kron and isopho- 
tal) magnitudes provide a good estimate of the total 
flux of an object (without the need for an aperture 
correction to be derived based upon the galaxy's as- 
sumed surface brightness distribution and the seeing, 
as is required to correct a small fixed aperture), the 
apertures must generally be large, containing areas of 
low surface brightness which increase the photomet- 
ric error. Also, an adaptive aperture might be a dif- 
ferent size in each measured band and thus sample a 
different component of the stellar population, poten- 

t webast.ast.obs-mip.fr/hyperz/ 



tially leading to unphysical colours. For these reasons, 
and because of the problems using adaptive apertures 
at faint magnitudes (where errors in determining the 
correct isophotal or Kron radii can exacerbate pho- 
tometric errors in low signal-to-noise images), we de- 
cided to use fixed aperture photometry for our SED 
fitting. 

However, fixed apertures have their own issues to 
consider. In the ideal situation, we would have a set of 
images with the same PSF, perfectly registered astro- 
metrically. It would then be possible to extract pho- 
tometry using apertures placed in exactly the same 
positions in each image, and we could be confident 
that, whatever the aperture size, we would be sam- 
pling the same fraction of the light from the same 
regions of each galaxy studied. But this is far from 
the case when combining data taken with different in- 
struments, under different seeing conditions, as we are 
doing here. It was undesirable on confusion grounds 
to degrade the optical and near-IR data to match the 
large PSF of the Spitzer bands, so we chose the aper- 
ture size to minimise the systematic differences in the 
fraction of the flux measured in each band. 

Clearly, there is a trade-off to be made — too 
large an aperture is vulnerable to confusion and in- 
creased background noise, while an aperture which 
is very small will contain very different fractions of 
the total flux for images with different PSFs. To in- 
vestigate this trade off in detail, we looked at the 
mean difference between aperture magnitudes and the 
SExtractor MAG_AUTO magnitudes, which uses an 
adaptive Kron aperture to estimate the total flux. For 
the smallest aperture available in all bands (2 arc- 
sec) there were large offsets (of order 0.2-0.3 mags) 
between the aperture and Kron magnitudes, but the 
most importantly these offsets were very different in 
the different bands. With the next available aperture 
diameter (3 arcsec) the mean differences in the aper- 
ture and Kron magnitudes were much smaller — less 
than a tenth of a magnitude. Deeming this an ac- 
ceptable systematic error, we decided to use 3-arcsec 
apertures for all photometry. 

To combine the catalogues across the optical and 
infrared bands, we assumed that the closest match 
between extracted objects, within a search radius of 
2 arcsec relative to the radio position (chosen to allow 
for differences in astrometry between different photo- 
metric bands), was the correct association, with the 
constraint that a catalogue entry was only permitted 
to match once per band. 

We used both empirical and evolving galaxy 
SEDs in the template fitting. The empirical templates 
used were the ones supplied with the HYPERz code, 
which are four typical observed galaxy spectra taken 
from Coleman, Wu & Weedman (1980; henceforth 
CWW), extended in the ultraviolet and infrared by 
matched models (Bruzual & Chariot 1993). 

The evolving templates were generated using 
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Galaxev (Bruzual & Chariot 2003; henceforth BC*), 
with 14 exponentially decaying star formation rates 
with e- folding timescales ranging from 0.1 to 30 Gyr, s 
plus an instantaneous burst, a finite burst and a n" 
constant star formation rate template, all with solar £ 
metallicity and a Salpeter initial mass function for 
the stars. This is a larger number of BC templates ~% 
than used by Bolzonella et al. (2000), but our inten- 
tion was to span the space of possible galaxy types as 
finely as possible, in the hope that there would always 8 
be a template with similar star formation history to o 
any object in the dataset, and that this would then o 
be the best fit. a 

Both empirical and evolving SEDs were used be- 
cause parallel work by van Breukelen et al. (2006) 
seemed to indicate that, while the evolving templates 
gave slightly more accurate results on spectroscopi- 
cally identified (and thus also reasonably bright) sam- 
ples, they more quickly broke down, giving catas- 
trophic redshift estimates, as signal-to-noise was re- 
duced. Catastrophically wrong redshift estimates pro- 
duced by the two sets of templates are generally differ- 
ent, so agreement between the results obtained with 
the CWW and BC SEDs is useful to improve confi- 
dence that a particular redshift estimate is correct. 
Where the estimates differ significantly, it is more of- 
ten the case that the CWW estimate is closer to real- 
ity, though this effect is difficult to quantify precisely 
since catastrophic redshift estimates are a function of 
the template set, as well as the precise signal-to-noise 
in each observed band. 

Figure Al presents the SEDs constructed for our 
radio sources; red colour represents the BC templates 
and black the CWW templates. We decided to use the 
BC derived value in the case where no spectroscopy 
is available. 

In Fig. 2 we compare 24 spectroscopic redshifts 
with the photometric redshifts (Table 2); optical spec- 
tra of these objects are presented in Fig. Al. Note that 
objects significantly deviating from the z spC c = %>hot 
line have unclear photometric redshifts (see Table 1). 
Three objects in our sample are classified as quasars 
(see Sec. 3.3), two of which fall into the z Bpcc > z phot 
category; the third, sxds_0009, may have a good pho- 
tometric redshift because the tentative broad-line 
emission could be scattered light, with the nucleus 
making an insignificant contribution at the longer 
wavelengths dominating the calculation of the pho- 
tometric redshift. 

The median redshift for our radio-source sam- 
ple, including photometric redshifts, is z = 1.095. 
The best fit line for the galaxies 'C gives an excel- 



t When running HYPERz with the CWW templates, 
slightly different apertures were used for the magnitudes. 
An exception is sxds_0016, where magnitudes were mea- 
sured directly from images and then the BC and CWW 
templates were fitted to the SED. 
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Figure 2. Comparison of spectroscopic and photometric 
(based on BC templates) redshift in the 24 (of 37) cases 
where both are available. Symbols indicate optical/near- 
infrarcd classification as adopted in Section 3.3, where 
filled symbols denote objects with spectroscopic redshifts: 
red filled circles for quasars 'Q', blue filled triangles for ob- 
scured quasars 'OQ', black filled squares for possible galax- 
ies 'G?', black squares for secure galaxies 'C, green upside- 
down triangles for starbursts 'SB', orange diamonds for 
weak quasars 'WQ' and light-blue stars for 'BL Lac' ob- 
jects (note that the 'BL' object is not present due to ab- 
sence of spectroscopic redshift). The light- blue dotted line 
corresponds to z spcc = ZphotBC- The solid black line fits 
the 15 galaxies 'C that have both spectroscopic and pho- 
tometric redshift: 2 spcc _G = 0.88 X z p hotBC-G — 0.03. 
Names in red represent objects with broad-line optical 
features. Outlying points are labelled with their source 
number. Objects with uncertain photometric redshifts are 
plotted with symbols, a factor 1.5 larger than the others. 

lent agreement between spectroscopic and photomet- 
ric redshifts, within the errors. The linear Pearson 
correlation coefficient is 0.89 and the Spearman cor- 
relation coefficient is p = 0.96 with 99.9% probability 
for a correlation. This makes us confident to use pho- 
tometric redshifts in the case where no spectroscopy 
is available. Consequently, we adopt the photometric 
redshifts in bold- type in Table 1. 

3.3 Optical/infrared radio-source 
classification 

Figure 3 shows the rest frame luminosity of the SXDS 
radio sources at 24 fim versus their redshift. The 
value of log 10 (Z/24Mm) is calculated using the spectral 
index measured from the 24 /jm data point to the 
nearest lower-A detection in the observed frame SED 
(Fig. Al). 

We use optical/IR observations (Fig. 3) to clas- 
sify a radio source as either Quasar 'Q', Obscured 
Quasar 'OQ', Galaxy? 'G?', Galaxy 'G', Starburst 
'SB', Weak Quasar 'WQ' or BL Lac 'BL', as described 
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Figure 3. Logarithm of the rest-frame Luminosity L24/jm 
at 24 /im versus redshift z. Symbols are the same as 
in Fig. 2. The dotted line corresponds to XL 24um = 
10~ 1,8 Z/Edd f° r a quasar with Mgjj ~ 10 8 Mq, as dis- 
cussed in Section 3.3. Outlying points are labelled with 
their source number. 

below. We deem that nuclear accretion is 'significant' 
in objects that obey log 10 (L24Mm/[WHz -1 sr -1 ]) > 
23.1 (or [AL] 2 4 M m > 10 37 ' 3 W)§. This value corre- 
sponds to [AL]24 M m > W~ 18 LEdd, a typical lower 
limit for quasars at z ~ 1 (see Fig. 2 McLure & 
Dunlop 2004), for a black hole mass Mbh > 10 s Mq, 
a typical lower limit for radio sources (McLure et 
al. 2004); L-^dd is the Eddington luminosity Lsdd = 
10 39 1 x (M/1O 8 M )W. 

We then define the following categories: 

i) Q: Broad lines in the optical spectrum (3/37 cases). 
None of these are detected at 24 /jm, although their 
limits are insufficient to rule out significant accretion. 

ii) OQ: Objects with a 24-pm detection (5/37 cases) 
and with sufficient L24 to represent significant accre- 
tion. In all cases we believe that the nucleus seen at 
24 fim is obscured optically because the spectra show 
no broad lines. This class may be incomplete in that 
some objects in the *G?' class, as described next, have 
limits above this critical value. 

iii) G?: A galaxy that has a 24-/xm limit consistent 
with it lying above the log 10 (Z/24fim/ [WHz _1 sr -1 ]) = 
23.1 line^ (11/37 cases). 

iv) G: All other objects (15/37 cases) without sig- 



3 Other authors use slightly different values to separate 
'luminous' quasars (and broad-line radio galaxies) from 
objects with lower accretion rates (e.g. Ogle et al. 2006 
adopt a value of [AL]i 5flm ~ 10 36 9 W). 
* Because of the 24 fim flux density limit these objects are 
the high-redshift (and hence, because of the 1.4-GHz flux 
density limit) high-Li .4GHz sub-set of the objects lacking 
Spitzer 24 (im detections. 



nificant accretion, unless they fall into three special 
categories defined by properties derived from consid- 
erations of the optical spectroscopy, the SED and the 
optical structure: SB: evidence from the SED and op- 
tical spectra of a starburst component (1/37 cases; see 
notes in Appendix A); WQ: evidence from the SED 
and optical spectra of a quasar component (see notes 
in Appendix A) but no 24 ^m detection (1/37 cases); 
BL: Featureless red continuum and a point source at 
K (1/37 cases). 

Objects labelled as 'Q', 'OQ' and 'G?' include all 
objects with, or potentially with, significant nuclear 
accretion, and are represented by filled symbols in 
Fig. 4. Note, however, that the definition of 'signifi- 
cant' accretion is somewhat arbitrary as illustrated by 
the case of the 'WQ' sxds_0025. Some accretion may 
be necessary in many, or perhaps all FRIs, and espe- 
cially those near the FRI/FRII divide, if the radio lu- 
minosity is assumed to be derived from the Blandford- 
Znajek mechanism (Cao & Rawlings 2004). 



4 DISCUSSION 

4.1 Accretion Indicators 

Inspection of the SEDs in Fig. Al shows that some 
of our objects have excess emission at 24 fj,m above 
that expected from extrapolation of the stellar pop- 
ulations. This is quantified via a measure of mid-IR 
excess, log 10 ([^L] 

lOptmrcst l\vL] 1 pi m rest 

) (Table 3). We 

chose a rest-frame value of 10 fj,m as our fiducial point 
as this corresponds to ~ 24 fim in the observed frame 
for redshifts close to the median redshift of our sam- 
ple. A comparison of mid-IR excess and blueness is 
presented in Fig. 4 where a positive correlation is 
clear. The generalised Spearman correlation calcu- 
lated using survival analysis statistics ASURV (Laval- 
ley et al. 1992) is 0.636 with a 99.9% probability for 
a correlation. 

Consider a simple model in which blueness is con- 
nected to mid-IR excess through the following equa- 
tions: 

= klX[^L]l (Jmres t+k2X[^L]l0amre a t =>(1) 
L J 4000Arcst l j ^ l j v / 

log 10 (e)x^x (gS==± ) +log 10 (k0», (2) 

Ki y [fJ-iJl/imrest J 

where ki encodes the contribution of the stellar pop- 
ulation of a passively evolving galaxy formed at high 
redshift (z > 5), and k2 the mid-IR-excess parame- 
ter that we are looking to calculate for this sample 
of radio sources. This model assumes that light from 
the nucleus with intrinsic optical luminosity L opt is i) 
absorbed by dust and re-emitted in the mid-IR gen- 
erating luminosity [vL]ioiimrest and ii) scattered, gen- 
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Figure 4. Bluoncss versus mid-IR excess for the 37 bright- 
est SXDS radio sources: Symbols indicate optical/IR clas- 
sification, as explained in Fig. 2 and Sec. 3.3. Table 3 
presents measurements of those quantities. The black dot- 
ted line corresponds to the best-fit line in the log-linear 
plane when all objects were treated as detections; the slope 
and intercept are 0.51 and -1.93 respectively, giving [see 
Eqn (2)] ki ~ 0.01 and k 2 ~ 0.01. The blue solid line 
corresponds to the best-fit line in log-linear plane where 
objects without detections at 24 fim were treated as limits; 
the slope and intercept are 0.09 and -0.34 respectively, giv- 
ing ki ~ 0.46 and k 2 ~ 0.10. The Buckley- James method 
in the survival analysis (ASURV) statistics package (Laval- 
ley et al. 1992) was used in these calculations. 'SB', 'BL', 
'Q' and 'WQ' objects were excluded from the fits as ex- 
plained in Sec. 4.1. 



erating a contribution to luminosity \uL] 

° L MOOOArest 

Objects labelled as 'SB', 'BL', 'Q' and 'WQ' have 
been excluded from the fit since they have SEDs domi- 
nated by different physical processes to those assumed 
in the model described by Eqns (1) and (2). Without 
spectropolarimetry it is difficult to determine whether 
any broad lines are seen as scattered light from the 
nucleus or they represent a direct view of the nu- 
cleus. The 24-/xm luminosity versus redshift plane of 
Fig. 3 shows that all three broad-line objects, as well 
as the 'OQ' objects, are consistent with being above 
or at the limit of L 24Mm > 10 23 1 WHz^sr" 1 , or 
[AZ,] 24 „ m > 10 37 - 3 W. 

Fig. 4 shows best-fit lines for two scenarios: 1) 
all objects were treated as detections (black dotted 
line), and 2) objects are treated as upper limits ac- 
cording to their 24 fim detection (blue solid line). We 
note that these 'fitted' lines are highly sensitive to 
the locations in the log 10 ([vL] 4 ^^^^/[uL]i limitm t) vs 



II Equation (2) relies on the ln(l + x) x approximation 
which is only accurate around and below the knees of the 
functions plotted in Fig. 4 
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Figure 5. Radio Luminosity at 1.4 GHz log 10 (Li.4GHz/ 
[WHz - 1 sr~ 1 ]) versus the radio spectral index a as in Table 
1: symbols are the same as in Fig. 2. The horizontal lines 
show the RLF and FRI/FRII breaks calculated from the 
Willott et al. (2003) and the Fanaroff & Riley (1974) values 
respectively using a typical steep-spectrum spectral index 
of 0.8. The vertical line at a = 0.5 shows the conventional 
border between steep and flat-spectrum sources. Objects 
that are discussed in the text (see Sec. 4.1) are labelled 
with their names. 
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lO^imrcst 



/\vlj\iiimrest) plane of a small num- 



ber of objects, and hence are not robust, although 
they certainly seem to delineate rough bounds to the 
data. Averaging these results we deduce ki ~ 0.2 and 
k 2 ~ 0.06, which agrees well with independent evi- 
dence for the properties of radio sources as we now 
argue. The value deduced for ki is in line with the 
expectations of template spectra of galaxies which 
formed their stars at high redshift. Optical polar- 
isation studies (e.g. Kishimoto et al. 2001) tell us 



that \uL] 

L J A 



0.01[^L] O pt, which is consistent 



400DArcst 

with our value of k 2 given that QSO SED studies 
suggest [^L]io M mrcst ~ 0.3[i/L] op t (Rowan- Robinson 
1995), and 0.06 x 0.3 ~ 0.01. We conclude that when- 
ever nuclear accretion is significant in our sample of 
radio sources, dust in the torus absorbs 30% of the 
photons and dust above and below the torus scatters 
>1% of the photons. 

In previous studies, the quasar fraction has been 
defined as the number of sources with quasar-like op- 
tical features (e.g. optical broad lines) and has a value 
of ~ 0.1 — > 0.4 over the relevant range of Li.4ghz (e.g. 
Willott et al. 2000). We introduce the term 'quasar- 
mode fraction' /qm to describe the fraction of objects 
with 'significant' accretion rates, as determined from 
the observed 24 /im flux density, to the total number 
of objects. 

In Figure 5 we present the radio luminosity at 
1.4 GHz of the SXDS radio sources against their ra- 
dio spectral index a (Table 1). A typical flat/steep- 



© 0000 RAS, MNRAS 000, 000-000 



14 Vardoulaki et al. 



spectrum radio spectral index a = 0.4/0.8, depending 
on whether its radio structure is compact/extended, 
was used in the case where the object was not ob- 
served at 610 MHz and is a limit at the 325-MHz ^ 
catalogue. Note that the median radio luminosity is £ 
Liaghz ~ 1.6 x lO^WHz - 1 sr -1 , which lies very close t 
to the FRI/FRII luminosity division. Crudely, the E 
SXDS sample is sensitive to objects above this di- 
vision at z>l and to objects below the FRI/FRII lu- . 
minosity break at z<l. 

We see from Fig. 5 that the single 'BL' object has 
a flat radio spectral index, as expected. It is clear that 
nearly all the other flat-spectrum object have signifi- 
cant accretion, although some exceptions do exist, e.g. 
object sxds_0028 lies close to the boundary of a = 0.5 
although it has significant errors in the calculation of 
its spectral index. Object sxds_0022 also lies close to 
the flat /steep a boundary, which is a typical value 
if star formation is making a significant contribution 
to the radio flux density (see notes in Appendix A). 
Object sxds_0025 is a 'WQ' with a flat radio spec- 
tra index and is expected to be in that region since 
we believe we are viewing the central region directly 
(see notes in Appendix A). The fact that it lies be- 
low the limit for 'significant' accretion is consistent 
with it being a weak quasar. Thus, the only clear ex- 
ception is sxds_0036, which is a 'G' and there is an 
uncertainty in the calculation of its spectral index. 
The only clear region where the quasar-mode frac- 
tion is very low, is for steep-spectrum objects below 
the FRI/FRII break, where / QM ~ 0.1. 

Figure 6** shows the radio luminosity at 1.4 GHz 
versus the projected linear size D. For SXDS sources 
that are compact we assume a limiting angular size 
9 — 2.0 arcsec (Simpson et al. 2006). We see that 
nearly all 'Q', 'OQ' and 'G?' objects of our sample lie 
above the 'FRI/FRII' luminosity break with the ex- 
ception of the 'OQ' sxds_0034 (the three 'G?' objects 
near sxds_0034 lie close to the boundary of significant 
accretion). Above the FRI/FRII break we find that 
/qm ~ 0.5 — 0.9 (the lower value assumes the 24 /im 
limits are much higher than the true 24 fj,m values, 
whereas the higher value assumes the true values lie 
just below the limits). The one clear exception in this 
regime is sxds_0001, which has no evidence of a QSO 
and a clear Twin- Jet (FRI) radio structure. 

The quasar-mode fraction drops dramatically be- 
low the FRI/FRII break, and whether or not one ex- 
cludes some of the compact (D < 100 kpc) sources as 
potentially part of a separate (beamed) population, 
then /qmS; 0.1 because nearly all objects are galax- 
ies 'G'. The clear counter-example here is sxds_0034, 
the only 'OQ' below the FRI/FRII break, and poten- 
tially an optically-obscured example of unobscured 



** In the similar figure in Vardoulaki et al. (2006) we used 
photometric redshifts to estimate values of £i.4GHz used 
in the plots. 



1 

: • Q 

: ▲ OQ 

: ■ G? A 

: □ G 
r 

: o WQ 

: * BL Lac • ▲ 

■ 


1 


- RLF Break 

■ ■ 


•at ' J 

A 


- FRI/II Break 'M 
7 ti34 
: V^D 

; % a a 
: *e 

1 


□ ■ w : 

CP : 
□ □ : 
□ : 

i 



1 10 100 1000 

D (kpc) 

Figure 6. Radio Luminosity at 1.4 GHz log 10 (Li.4GHz/ 
[WHz - 1 sr -1 ]) versus largest projected linear size D. Sym- 
bols are the same as in Fig. 2. The horizontal lines show 
the RLF and FRI/FRII breaks calculated from the Willott 
et al. (2003) and the Fanaroff & Riley (1974) values respec- 
tively using a typical steep-spectrum spectral index of 0.8. 
Objects that are discussed in the text (see Sec. 4.1) are 
labelled with their names. 



FRI QSOs already studied in this radio luminosity 
regime (e.g. Sarazin et al. 1999, Blundell & Rawl- 
ings 2001). This object is spectroscopically confirmed 
at redshift z — 1.095. An excess of emission at 24 
/jm places it securely in the 'OQ' classification. This 
object is not a spectroscopically confirmed quasar 
but has high-excitation narrow lines in its optical 
spectrum and is compact both optically and at K; 
its radio map shows a seemingly unresolved source. 
This seems to be an optically-obscured 'radio-quiet' 
quasar in the radio luminosity regime in which most 
radio sources are thought to have FRI radio struc- 
tures. Unobscured 'radio-quiet' quasars are known in 
this regime ( see Heywood et al. 2007 and references 
therein), and claims have been made for obscured 
examples (Martinez-Sansigre et al. 2006). However, 
as found by the study of Martinez-Sansigre et al. 
(2006), the surface brightness sensitivity of existing 
radio data is insufficient to map any extended struc- 
ture at high redshifts. We also note the presence of 
broad absorption lines in 2 out of our 3 'Q' category 
objects. This is a high fraction compared to general 
QSO samples (e.g. Richards et al. 2003 find this frac- 
tion to be between ~ 0.03 and ~ 0.2 for optically 
selected samples of quasars). 

Both of Figs 5 & 6 suggest that we have a larger 
fraction of flat, compact radio sources in our sam- 
ple since it is selected at higher frequency than low- 
frequency selected radio samples like the 3CRR, 6CE 
and 7CRS: 51% ± 7% objects of our sample have 
D < 25 kpc, whereas this fraction is ~ 25% in 7CRS 
(e.g. Fig. 18 of Blundell, Rawlings & Willott 1999). 
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The quasar-mode fraction is generally high for fiat- 
spectrum objects independent of radio luminosity. 
Steep-spectrum and extended, presumably FRI-type 
radio sources (i.e. objects below the FRI/II break at 
a > 0.5) show no signs of significant accretion with 
just one exception. Cao & Rawlings (2004) studied a 
sample of 3CR FRI radio sources and found that the 
L/L E dd ratio in FRI-type radio sources is more than 
two orders of magnitude less than our limit for signif- 
icant accretion, which is in line with our findings. 

We have discovered that above the FRI/FRII 
division the quasar-mode fraction is of the order 
/qm ~ 1, where the median redshift of these objects 
is z m ed ~ 1.6. Below the FRI/FRII break /qm — > for 
extended steep-spectrum sources, where 2 mc d ~ 0.65. 
A key question is whether the presence of significant 
accretion is linked to the radio luminosity L ra d (or ra- 
dio structure) or to the redshift z. We note that from 
studies of the 3CR sample (in which the FRI/FRII 
division is at redshift z ~ 0.1 rather that z ~ 1, for 
our sample) that a similar division of mid-IR proper- 
ties occurs at the same FRI/FRII luminosity division 
(Ogle et al. 2006). We conclude that it is a luminos- 
ity (or structural) effect rather than an epoch effect, 
and in this paper we have shown that such an ef- 
fect persists from z ~ 0.1 to z ~ 1. We will attempt 
to de-couple the effects of radio luminosity and radio 
structure in a future paper. 



5 CONCLUSIONS 

We have presented a study of the brightest SXDS 
radio sources with flux density at 1.4 GHz above 2 
mjy. Our sample is not as yet spectroscopically com- 
plete (24 out of 37 have spectroscopic redshifts), but 
Fig. 2 gives us the confidence in the use of photomet- 
ric redshifts calculated with HYPERz using Bruzual- 
Charlot templates. The median redshift of the entire 
sample is z me( j ~ 1.1, which from Fig. 5 corresponds 
to objects close to the FRI/FRII division. Optical and 
infrared imaging, as well as radio observations are in 
place to give an insight to the type of radio sources 
at the bright flux density end of the 1.4-GHz SXDS 
radio survey of the SXDF. 

We have classified our radio sources as either 
quasars 'Q', obscured quasars 'OQ', possible galax- 
ies 'G?', secure galaxies 'G', starbursts 'SB', weak 
quasars WQ' or BL Lac objects 'BL'. Objects 'Q', 
'OQ' and 'G?' have significant nuclear accretion, 
since they obey log 10 (L24/am/[WHz _1 sr -1 ]) > 23.1 
(or [AL] > 10 37 ' 3 W). This value corresponds 
to [AL] 24nm > 10 1 ' fl LEdd for a black hole mass 
M BH > 10 8 M Q . 

The 37 brightest radio sources in the SXDS ex- 
hibit the following properties: 

• Flat-spectrum, presumably mostly Doppler - 
boosted, radio sources show evidence of accretion in 
almost all cases, independent of radio luminosity. 



• Significant accretion is also nearly ubiquitous above 
the 'FRI/FRII' luminosity break, i.e. at redshift z>l 
for objects in our sample. The quasar-mode fraction, 
defined here as the fraction of objects with significant 
accretion rates to the total number of objects - and 
not as the number of sources with quasar-like optical 
spectra - is ~ 0.5 — > 0.9. There are counter exam- 
ples like the low-accretion FRI sxds_0001, but they 
are rare. 

• At z<l most of the 'FRI'-regime sources have low 
accretion rates and /qmS; 0.1, but one high-accretion- 
rate object (sxds_0034) does exist; this is an obscured 
'radio-quiet quasar' which may have FRI-like jets. 

• The presence of significant accretion is linked to the 
radio luminosity L ra( j (or radio structure) of the radio 
source at z ~ 1, just as it is at z ~ 0.1. 

• Mid-IR excess and blueness are correlated, where 
the Spearman's generalised p = 0.636 with a 99.9% 
probability for a correlation. 

• This correlation is consistent with a model in which 
~ 30% of the QSO light is absorbed by a dust torus 
and >1% is scattered by dust above and below the 
torus. 

Study of the brightest SXDS radio sources in the 
mid-IR has revealed that the ratio of objects with sig- 
nificant accretion rates to the total number of objects 
is high for flat-spectrum radio sources, and for steep- 
spectrum sources above the FRI/FRII radio luminos- 
ity divide. This ratio is much lower for objects below 
the luminosity corresponding to the FRI/FRII divide. 
Future work on the brightest SXDS radio sources is 
likely to involve studies of radio structural proper- 
ties, host galaxies and other indicators of accretion 
like emission lines and X-rays. 
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APPENDIX A: NOTES ON THE 37 
BRIGHTEST SXDS RADIO SOURCES 

sxds_0001: The brightest radio source of our sam- 
ple is a twin-jet (TJ) galaxy at redshift z = 0.627. 
Spectroscopy was obtained both with VIMOS on the 
VLT and ISIS on the WHT. In the latter case the 
slit width was 2.5 arcsec and the position angle PA 
= 320°. This object is resolved at K. The NVSS flux 
density is higher than the SXDS value (Table 1), in- 
dicating the existence of diffuse radio structure sur- 
rounding the radio source. This is a typical galaxy 
'G', according to Sec. 3.3, with a prominent [Oil] line 
in its optical spectrum and no hint of 24-/im emis- 
sion. The apparent flattening between 610 MHz and 
325 MHz maybe due to inaccurate flux density entry 
in the 325-MHz catalogue. 

sxds_0002: This object has a classical double 
radio structure, as seen in the A-array radio map. 
It is also an obscured quasar 'OQ', since it has 
log 10 (Z/24Mm/W) > 37.3 and is resolved, and pre- 
sumably dominated by starlight at K. Optical spec- 
troscopy was obtained at the WHT using a slit with 
a width of 2.5 arcsec and a PA = 305°. The BRz- 
radio overlay in Fig.Al shows the possibility of merg- 
ing activity between the radio-source host galaxy and 
a nearby object. 

sxds_0003: This is a classical double radio 
source since the A-array radio map reveals com- 
pact FRII-like hotspots associated with an obscured 
quasar at redshift 2.433. Spectroscopy was obtained 
both with VIMOS at the VLT and ISIS at the WHT, 
where in the latter case, the slit width was 2.5 arc- 
sec and the PA = 313°. The optical spectrum lacks 
broad emission lines, a result that confirms our classi- 
fication of this radio sources as 'OQ'. As in sxds_0002, 
the BRz-radio overlay reveals possible galaxy-galaxy 
merging activity and it has a 24- /im excess in its SED. 
The Lya line in the emission spectrum of this object 
is spatially extended over ~ 3.4 arcsec ~ 25 kpc. 

sxds_0004: A compact flat-spectrum radio 
source with a SED rising to 24 fim. According to our 
classification scheme of Sec. 3.3, this is a BL Lac ob- 
ject, which is consistent with it having a flat radio 
spectral index (a = 0.164). The difference between 
SXDS and NVSS flux densities (Table 1) indicates 
that the radio-source flux density might be time vari- 
able. This object was observed using AAOMEGA at 
the AAT, but we only have a 1-D spectrum available, 
with one highly questionable emission line at 5900 A, 
and this does not give a reliable spectroscopic red- 
shift; we prefer thus the photometric redshift. 

sxds_0005: A compact steep-spectrum radio 
source and a galaxy without a secure spectroscopic 
redshift, since observations with the WHT didn't give 
a sufficient S/N to measure an accurate redshift (slit 
width 2.5 arcsec and PA = 310°). The optical spec- 
trum shows a very faint continuum, as expected from 
its photometry (Simpson et al. 2006) and the fact that 



the photometric redshift, z — 1.610, is in the 'redshift 
desert' where there is a lack of any strong emission 
lines. 

sxds_0006: A compact radio source that is pro- 
jected close to a bright star and hence without any 
entry in the mid-IR Spitzer catalogue. This object 
has no optical spectroscopy, but we are confident that 
the photometric redshift is accurate due to a good fit 
of the template to the SED (Fig. Al). It has an in- 
verted spectral index, as calculated from the 1.4 GHz 
and 610 MHz data; the limit at 325 MHz verifies this 
result. 

sxds_0007: A large (123 arcsec) FRII classical 
double radio galaxy at z = 1.325. The NVSS flux 
density is the sum of the integrated fluxes of the two 
lobes, which is the same value as the one measured 
from the NVSS map itself. Both 610 MHz and 325 
MHz data suggest a very steep radio spectral index 
for this source. 

sxds_0008: A compact radio source in a z — 
0.458 blue galaxy. According to our classification of 
radio sources, this is a starburst 'SB' with a star for- 
mation rate SFRfir ~ 100 Moyr -1 , calculated from 
the SED, using a Salpeter IMF and the star formation 
rate (SFR) relations of Kennicutt (1998). We have a 
good agreement of spectroscopic and photometric red- 
shifts. From the BRz-radio overlay of Fig. Al we see 
obvious galaxy-galaxy merging activity and a 24-^im 
excess in its SED. We believe the 24 /xm excess is more 
readily attributable to a starburst than an AGN be- 
cause the optical spectrum has bright [Oil] emission 
but no corresponding high-excitation lines. 

sxds_0009: A compact radio source missing J— 
and K — band data since it lies outside the UDS cov- 
erage (Fig. 1). The weak extensions is the B-array 
radio map (Fig. Al) to the east and the south-west 
are probably artifacts. The high NVSS flux density 
is due to confusion with a nearby source. This object 
is probably at z = 3.570, which agrees with the pho- 
tometric redshift (Table 1). The photometry of this 

object in the blue, B = 26.3 mag^, indicates that 
the apparent blue tail in the optical spectrum is a 
spurious feature. This object has a flat radio spec- 
trum and some evidence that a significant fraction of 
the optical light comes from a point source although 
there are no broad lines or a 24-/im excess. We classify 
it as a 'Q?' because of hints of broad wings to the Lya 
line and the significant blue continuum, although the 
limit at 24-/im is insufficient to confirm a QSO-likc 
mid-IR excess. 

sxds_0010: A compact radio source and a galaxy 
'G'. The HYPERz photometric code gives a good fit 
to the SED when we use the BC templates. There are 
no spectroscopic data available for this source. 

sxds_0011: A radio galaxy at z = 0.647 with TJ 
structure. The angular size was measured from the ra- 

tt (www.astro.livjm.ac.uk/ ~ cjs/SXDS/radio/) 
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dio overlay from the core to the largest extent of the 
north-east jet and separately to the south-east radio 
components, with the two values added together. The 
610 MHz flux density was measured around the same 
area as the 1.4 GHz flux density giving a radio spec- 
tral index between these values of 0.62. The NVSS 
flux density, though, suggests that there is extra emis- 
sion at large angular sizes. We adopt a typical steep 
radio spectral index of a = 0.8 since steep spectral 
indices are typical for objects with a FRI radio struc- 
ture; the spectral index calculated from the SXDS 
and 325 MHz data suggests an inverted spectral in- 
dex which is inconsistent with the fact that this is a 
large twin-jet radio source. We suggest that the 325 
MHz data are inaccurate in this case. 

sxds_0012: A large FRII galaxy at 2 = 0.865. 
The 325-MHz catalogue entry includes only part of 
the source, giving a nonsensical radio spectral index. 
We use the total 610-MHz flux density from the radio 
map to calculate the spectral index, whose steep value 
agrees with the CD radio structure. There are no signs 
of a QSO nucleus. 

sxds_0013:A spectroscopically confirmed quasar 
at z — 2.713 and a point source in the optical. The 
optical spectrum shows strong absorption features as- 
sociated with the broad lines. Such broad-absorption- 
line (BAL) QSOs are rare, but not unknown, in radio- 
selected samples. The B-array radio image shows a 
compact source. From the optical overlay of Fig. Al it 
looks like a galaxy is dominant in the K— band, since 
it is not a point source at K, although it is one in the 
optical. The flat radio spectral index, a — 0.08, agrees 
with the fact that this is a compact radio source. 

sxds_0014: This object appears to be red with 
a rapid drop in its SED from red to blue. In the 
if— band image, although it has a low S/N ratio, 
the source looks unresolved indicating that the radio 
source might be compact or there is a hidden quasar 
nucleus. The limit at 24 /»m is insufficient to rule out 
that this might be an obscured quasar 'OQ'; we define 
it as a galaxy 'G?' (see Sec 3.3). In Vardoulaki, Rawl- 
ings & Simpson (2006) we used 36 radio sources from 
this sample, accidentally excluding this one; all the 
subsequent identification numbers have moved down- 
wards in order to include this radio galaxy. The in- 
verted radio spectral index is consistent with the fact 
that this is object has compact radio structure. The 
tentative emission feature around A ~ 6000A in the 
optical spectrum of this object cannot provide us with 
a reliable spectroscopic redshift. 

sxds_0015: A possible 'double-double' radio 
galaxy. Such sources (Schoenmakers et al. 1999) are 
thought to be jets that have re-started, sometimes 
(as perhaps in this case) with a change of jet orien- 
tation, perhaps due to merging of two supermassive 
binary black holes (Liu et al. 2003). From the BRz- 
radio overlay of Fig. Al we see evidence of galaxy- 
galaxy merging in the form of tidal tails (Johnston et 
al. 2004), and there is also excess emission at 24 fim 



as can be seen in its SED. We use the radio spectral 
index calculated between 610 MHz and 1.4 GHz data, 
which gives a steep value that agrees with the FD ra- 
dio structure. We don't trust the 325 MHz data since 
they give a spectral index inconsistent with the radio 
structure. Although this object is defined as 'G', its 
detection at 24 /xm might indicate it is either a 'SB' 
or a 'WQ'. 

sxds_0016: A classical-double radio source, pos- 
sibly at z = 1.345. Shorter wavelengths in the spec- 
trum are, due to poor resolution, affected by the over- 
lapping spectrum of a nearby object. Photometry of 
this object was directly measured from images using 
sxds_0034 as a boot-strap calibrator, since it was not 
matched with any source within a 2" search radius. 
Using these data we give the photometric redshift in 
Table 1 with an associated uncertainty. According to 
these measurements, the object has excess emission 
at 24 pirn, which classifies it a 'OQ' according to Sec- 
tion 3.3 and Fig. 3. We use the spectral index calcu- 
lated between the 610 MHz and 1.4 GHz flux densi- 
ties, which agrees with the CD radio structure of this 
object. 

sxds_0017:A spectroscopically-confirmed broad- 
absorption- line quasar at z — 2.583, where the Lya 
line is heavily absorbed. The spectrum appears pecu- 
liar, since the CIV and CIII] lines have narrow cores 
and a red and blue broad wing respectively. The ra- 
dio map reveals a prominent hotspot to the north- 
east; the component to the south-west of the quasar 
is probably the weak side of an asymmetric double. 
We use the 325 MHz flux density to calculate a spec- 
tral index since it was not observed at 610 MHz. 

sxds_0018: A large FD (no compact hotspots in 
the A-array map) radio galaxy at z = 0.919. The ra- 
dio spectral index is calculated between the 610 MHz 
and 1.4 GHz data, where the steep value is consistent 
with the extended FD radio structure. The 325-MHz 
catalogue entry is probably inaccurate, since plausi- 
ble extrapolation of the 610 MHz flux density suggests 
that the object has more flux density at 325 MHz than 
the limit of that catalogue. 

sxds_0019: A FD (no compact hotspots in the 
A-array map) radio galaxy at z = 0.694. It has a 
very asymmetric radio structure. We assume a = 0.8 
due to the extended FD radio structure; the inverted 
spectral index calculated between the 325-MHz and 
1.4-GHz data is inconsistent with the radio structure 
of this object. 

sxds_0020: This object has an unusual radio 
structure that we classify, uncertainly, as FD. The 
steep spectral index calculated between 610 MHz and 
1.4 GHz is consistent with the extended radio struc- 
ture. The 325-MHz entry is inaccurate since both the 
radio structure and the 610-MHz to 1.4-GHz spectral 
index indicate that it has a steep radio source. There 
are no spectroscopic data available for this object and 
the photometric estimate is highly uncertain. 



© 0000 RAS, MNRAS 000, 000-000 



The 37 brightest radio sources in the SXDF 19 



sxds_0021: A radio galaxy 'G' and a compact 
radio source. The radio structure to the north-west 
side of the object is probably not associated with the 
brightest radio source as it has a separate galaxy iden- 
tification. The photometric redshift is used since there 
is no spectroscopy available. 

sxds_0022: A compact radio source associated 
with a galaxy at z = 0.248; note the spectacular blue 
arc which is probably a remnant of star formation as- 
sociated with a galaxy-galaxy interaction. The south- 
ern and eastern radio components are probably ar- 
tifacts. This is possibly a post-starburst galaxy, or 
an 'E+A', i.e. an elliptical with additional absorption 
lines (e.g. H/3 in Fig.Al) in their optical spectra due 
to stars seen ~ 10 8 yr after the starburst has finished 
(Goto 2007). Star formation might be contributing 
to the radio emission. There is some evidence for ex- 
cess emission at 24 /jm as well as dramatic tidal tails 
(Fig. 1). We classify this radio source as a 'G'. 

sxds_0023: A compact radio galaxy at z = 
0.586. The A- array map reveals a possible double 
structure, although it might be the base of a FRI- 
like jet. There is a fairly good agreement between the 
spectroscopic and photometric redshifts. 

sxds_0024: A TJ radio galaxy at z = 0.516. The 
radio spectral index that is calculated between 325 
MHz and 1.4 GHz is inconsistent with the TJ radio 
structure of the source; so we use 610 MHz flux den- 
sity. 

sxds_0025: A compact weak quasar 'WQ' prob- 
ably at z — 0.9636 if the single marginally broad line 
in the spectrum is Mgll (the deconvolved FWHM is ~ 
900 kms -1 ). We are not certain that the [OII]A2470A 
emission line is real, since spectroscopic studies show 
that this is quite a weak feature (e.g. Serjeant et al. 
1998). The photometric redshift is inaccurate due to 
the absence of J— and K— band data (Fig. 1) and the 
fact that the optical continuum appears to be dom- 
inated by blue QSO light; the rest-frame EW of the 
putative Mgll line is ~ 50 A and thus within the 
spread seen in directly viewed quasar nuclei (Francis 
et al. 1999). The templates used in the HYPERz code 
cannot provide a proper fit to the SED. The NVSS 
flux density is higher than expected due to confusion 
by sxds_0027. This object was not observed at 610 
MHz since it is outside the coverage of the survey. 
We assign a typical flat radio spectral index a — 0.4 
consistent with the compact radio structure; we do 
not trust the 325-MHz catalogue to be complete. 

sxds_0026: A TJ radio galaxy at z = 0.450. The 
companion galaxy to the west is probably not produc- 
ing radio emission. According to the radio spectral in- 
dex calculated between 610 MHz and 1.4 GHz data, 
this is a steep-spectrum radio source. The fact that 
the 325-MHz entry is a limit suggest that we should 
not trust that the 325-MHz catalogue is complete. 

sxds_0027: We classify this object as 'OQ' since 
it has excess emission at 24 pirn. We use the photo- 
metric redshift since the optical spectrum shows a fea- 



tureless continuum. The longer-wavelength SED looks 
more like a BL Lac since it can easily be fitted by 
a power-law with very steep spectral index a ~ 2, 
and is consistent with the flat continuum in the op- 
tical spectrum of this object. The 'tails' in the radio 
map are probably artifacts, so its compact nature is 
consistent either with it being a BL Lac object or 
an 'OQ'. The object is resolved at K which is why 
we prefer the 'OQ' category. This radio source is not 
observed at 610 MHz. We assign a typical flat radio 
spectral index a — 0.4 consistent with the compact ra- 
dio structure; we do not trust the 325-MHz catalogue 
to be complete at the flux density level appropriate 
to this source. The NVSS flux density is higher than 
expected since with the NVSS resolution we can only 
marginally distinguish this object from the nearby ob- 
ject sxds_0025. 

sxds_0028: A compact radio galaxy, probably 
at z = 0.627. The radio map shows a seemingly dou- 
ble source, but in reality there are two radio sources, 
as there are good optical IDs under both apparent 
radio lobes. The object under the north-west 'lobe' 
is sxds_0028, whereas the one under the south-east 
'lobe' is sxds_0042 (a compact radio source at redshift 
z — 0.382, which is just below our flux density limit). 
The flattening of the radio spectral index between 
610 MHz and 325 MHz suggests that the 325-MHz 
catalogue is incomplete; the 325 MHz flux should be 
higher than the limit of that catalogue. 

sxds_0029 A compact radio source that looks 
unresolved in the K— band. Due to absence of spec- 
troscopic data, we use the photometric redshift, which 
suggests a fairly good fit to the SED (Fig. Al). 

sxds_0030: This probably is a TJ radio galaxy 
at redshift z = 0.535 with the jets swept back forming 
a narrow-angle tail. There is a fairly good agreement 
between the photometric and spectroscopic redshifts. 

sxds_0031: This is a compact radio galaxy. The 
radio component to the west is probably not associ- 
ated with this source as it has its own separate galaxy 
identification. The fairly good fit of the photometric 
redshift makes us confident to use it in the absence of 
spectroscopic data. 

sxds_0032: This is probably a FD radio galaxy; 
no A-array data available to confirm lack of compact 
hotspots. Without any spectroscopic redshift, the BC 
templates provide a very good fit to the photometric 
data and a believable photometric redshift. The 325- 
MHz catalogue misses this object. According to the 
spectral index calculated between the 610 MHz and 
1.4 GHz data, this object should have a flux density 
at 325 MHz that would give a steep spectral index. 

sxds_0033: A TJ radio source at z = 0.647. It 
belongs to a small group of galaxies that is interacting 
and probably merging with a larger group of galaxies 
(Geach et al. 2007). The feature at 5426 A in the 
optical spectrum (Fig. Al) is CIII] emission from a 
lensed background object at z = 1.847 (see Geach et 
al. 2007). It is not clear whether the excess emission 
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at 24 fim is associated with the lens or the galaxy. 
However, the majority of the radio emission probably 
originates in the lcnsing source. The steep spectral 
index calculated between 610 MHz and 1.4 GHz data 
indicates that the 325 MHz flux density should be 
higher than the catalogue entry. 

sxds_0034: A compact radio source at z = 1.095 
with a clear QSO component at the wavelengths 
probed by Spitzer data, a point-like optical structure, 
and high-excitation narrow lines in its spectrum. It 
is, most plausibly, a Doppler-boosted object with low 
extinction and possible time variability. This is the 
one object with significant accretion that lies below 
the FRI/FRII break (Fig. 6), indicating that it is a 
FRI 'radio-quiet' quasar. 

sxds_0035: A compact radio source and a galaxy 
without spectroscopy or J— /K— band data. The HY- 
PERz code fits the 4000 A break between J and K, 
which gives a photometric redshift of z = 2.715. 

sxds_0036: A compact radio source and a galaxy 
'G' at z = 0.871. The 610 MHz flux density is an 
upper limit, giving an uncertain measurement for the 
radio spectral index. 

sxds_0037: A compact radio source and a galaxy 
'G' that is not observed in the J and K band (Fig. 1), 
and which therefore has a tentative photometric red- 
shift. This radio source was outside the coverage of 
the 610 MHz survey. We use a typical flat radio spec- 
tral index a = 0.4, consistent with the compact radio 
structure of this source; we don't trust the 325-MHz 
catalogue to be complete at this flux density level. 
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Figure Al. Top: if -band images and overlaid radio con- 
tours from the VLA for sxds_0001; either the B/C-array 
(beam ~ 4 arcsec 2 ) or the A/B/C-array map (beam ~ 1.8 
arcsec 2 ) is shown, with the image labelled as B or A re- 
spectively. For sxds.0009, sxds.0017, sxds_0025, sxds_0027 
and sxds_0037 there are no A"-band images since these ob- 
jects are beyond the edge of the UDS survey (see Fig. 1), so 
we use z-band images instead. The radio structural classi- 
fication is shown on the bottom right of the image. Second 
from top: Overlays of B— , R— and 2— band images and 
B-array radio maps from Simpson et al. (2006). Second 
from bottom: Observed frame SEDs of our objects, where 
we also give the best-fit SEDs, their spectroscopic red- 
shift z sp cc (where applicable) and their photometric red- 
shift calculated from HYPERz: the red SED is calculated 
using BC templates and the black with CWW templates 
(see Section 3.2). In all the overlays, North is up and East 
is to the left. The contour levels follow: 2 ncontour_1 X a, 
where n con tom-s is the number of contours and o the noise 
level. Bottom: Optical spectra of 25 of the 37 SXDS radio 
sources. The character A' denotes atmospheric absorption 
feature, 'S' a featured caused by the dichroic, 'C a cosmic 
ray, and 'F' fringing on the CCD. All spectra have been 
smoothed over 3 pixels. The following spectra are not flux 
calibrated: sxds_0009, sxds_0028, sxds_0030. 
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Figure Al. The rest of Fig. Al can be found at the URL: www — astro.physics.ox.ac.uk/ ~ eleniv / figure_Al.html 
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